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Abstract. Measurements of the inclusive yields of n° mesons in p+p and d+Au collisions at center 
of mass energy ^/snn — 200 GeV and pseudorapidity (tj) = 4.00 (d beam direction) are reported. 
The yield for p+p collisions is in general agreement with perturbative QCD calculations. The d+Au 
yield is in agreement with a calculation which models the Au nucleus as a Color Glass Condensate 
for forward particle production. The nuclear modification factor derived from the inclusive yields is 
qualitatively consistent with models which suppress the gluon density in nuclei. 
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In contrast to the nucleon, very little is known about the density of gluons in nuclei 
HI]. For protons, the gluon parton distribution function (PDF) is constrained primarily 
by scaling violations in deeply-inelastic lepton scattering (DIS) measured at the HERA 
collider [2] . The data are accurately described by QCD evolution equations that allow 
the determination of the gluon PDF As the momentum fraction of the parton (xsj) 
decreases, the gluon PDF is found to increase. At sufficiently low xgj, the increase in 
gluon splitting is expected to become balanced by gluon-gluon recombination, resulting 
<D ! in gluon saturation. In nuclei, the density of gluons per unit transverse area is expected 

to be larger than in nucleons, and so is a natural environment in which to establish if, 
and under which conditions, gluon saturation occurs. Quantifying if gluon saturation 
occurs at RHIC energies is important because most of the matter created in heavy ion 
collisions is expected to evolve from an initial state produced by the collisions of low- 
^ ■ xgj gluon fields in the nuclei |2D. Fixed target nuclear DIS experiments are restricted in 

the kinematics available and have determined the nuclear gluon PDF for xgj ~ 0.02 JH]. 
In a conventional pQCD description of d(p)+A collisions at Jsj^ = 200 GeV, inclusive 
forward hadroproduction probes the nuclear gluon density over a broad distribution of 
xgj peaked around 0.02 and extending well below xgj ~ 0.005 0]. 

Using factorization in a perturbative QCD (pQCD) framework, the PDF's and frag- 
mentation functions (FF's) measured in electromagnetic processes can be used in the de- 
scription of hadronic scattering processes. In p+p collisions at y/s = 200 GeV, factorized 
leading twist pQCD calculations have been shown to quantitatively describe inclusive 
n° production over a broad rapidity window . In pQCD, forward % production in 
p+p collisions probes low-xg ; gluons (g) in one proton using the valence quarks (q) of 
the other. Recently, the yield of forward negatively charged hadrons (h~) in the d-beam 
direction of d+Au collisions was found to be reduced when normalized to p+p collisions 
0]. The reduction is especially significant since isospin effects are expected to suppress 
h~ production in p+p collisions, but not in d+Au collisions flj]. 

Many models attempt to describe forward hadroproduction from nuclear targets. 
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FIGURE 1. Inclusive n production cross section for p+p collisions. The error bars, often smaller 
than the points, combine the statistical and point-to-point systematic errors. (Left) Versus pion energy 
(E K ) at fixed pseudorapidity (77). The curves are NLO pQCD calculations using two sets of fragmentation 
functions (FF). (Right) Versus transverse momentum (pj) at fixed Feynman-x (xp). Curves here and in 
the inset are fits to the data of the form shown on the plots. (Inset) Versus (1 — xp ) at pp — 2GeV/c. 



Saturation models 18] include a QCD based theory called the Color Glass Condensate 
(CGC) Ja, llOll . Another approach models quarks scattering coherently from multiple 
nucleons, leading to an effective shift in the xbj probed [11]. Shadowing models suppress 
the nuclear gluon PDF in a standard factorization framework lfl2ll . Parton recombination 
models modify the fragmentation of a quark passing through many gluons II 1 3fl . Other 
descriptions include factorization breaking in heavy nuclei II14I1 . More data are needed 
to constrain the mechanisms for forward hadroproduction in nuclear collisions. 

In this paper we present the inclusive yields of high energy n° mesons at (rj) = 4.00 
(77 = — ln[tan(0/2)]) in p+p and d+Au collisions at ^/smn = 200GeV. The results 
are compared with theoretical predictions. The rj dependence of the normalized cross 
section ratio is also presented. Analysis of the azimuthal correlations of the forward n° 
with coincident hadrons at midrapidity is presented elsewhere H15I1 . 

A forward n° detector (FPD) was installed at the Solenoidal Tracker At RHIC (STAR) 
to detect high energy n° mesons with 3.3 < rj < 4.1. In the 2002 run, p+p collisions at 
y/s = 200 GeV were studied with a prototype FPD o]. In the 2003 run, p+p collisions 
were studied with the complete FPD and exploratory measurements were performed 
for d+Au collisions at ^/snn = 200 GeV. Details of the FPD performance can be found 
in Ref. Ill 611 . The luminosity was determined by measuring the transverse size of the 
colliding beams and the number of colliding ions, giving a normalization error of ~ 1 1%. 

The cross sections for p + p — » n° + X at (77) = 3.3, 3.8, and 4.00 are presented in 
Fig. □(left) |H [3]. The ^-dependent systematic error at (7]) = 4.00 is 10-20%, 
dominated by the energy calibration of the FPD. The absolute 7] uncertainty contributes 
11% to the normalization error iflfjll . The curves are NLO pQCD calculations using 
the KKP and Kretzer sets of FF's differing primarily in the gluon-to-pion FF, Dg. At 
7] = 3.3 and 3.8, the data are consistent with KKP. As pj decreases further, the data 



begin to undershoot KKP and approach consistency with Kretzer. This is consistent with 
the trend seen at midrapidity Jg]. At low pp, the dominant contribution to 7T° production 



becomes gg scattering, making D n g the dominant FF II 1711 . In Fig. [T] (right), the cross 
section is shown versus pp at fixed Feynman-jc (xp 2E n /\/s), and versus (1 — xp) at 
px = 2 GeV/c (inset). As was reported at the CERN ISR 111 8ll . the forward rapidity yields 
are rapidly changing functions of both xp and pp. 

The cross section for d + Au -> 7t° + X at ^/s^ = 200 GeV and (rj) = 4.00 is 
presented in Fig. El (left). The .E^ -dependent systematic error is « 20%, dominated by 
the background correction because, on average, 0.5 more photons are observed in d+Au 
than in p+p collisions per event with > 30 GeV detected in the FPD. The curves are 
LO calculations from Ref. II Oil , using CTEQ5 PDF's and KKP FF's convoluted with 
a dipole-nucleus cross section which models parton scattering from a CGC in the Au 
nucleus. The "MV" and "No DGLAP" curves neglect QCD evolution of the Au wave 
function and the PDF/FF, respectively. The difference in the slopes with and without 
evolution is greater than the slope change from LO to NLO for p + p — ► n° + X at 
r\ = 3.8. The calculations in Fig.EJare normalized by a pp -independent ^-factor of 0.8. 
This is smaller than the ^-factor used to normalize the theory to the d + Au — > h~ yield 
at the nominal value of T] = 3.2, in the same direction as the renormalization needed to 
scale the NLO calculations with KKP FF to p + p — > 7t° data at rj = 4.00. Note from 
Fig.[T]that a change of Arj « 0.05 at these kinematics results in Aa/a ~ 35% at fixed 
xp. The slope of the n° yield is consistent with the calculation where the PDF and FF 
evolve a la DGLAP, and includes small-x^/ evolution of the Au wave function. 

Nuclear effects on particle production are quantified by the nuclear modification 
factor, R^ Au , which can be defined for minimum-bias events by the ratio, 



>Y _ Ed 3 o/dp 3 (d+Au^Y+X) 
<AW<dr Ed*o/dp*(p + p-+Y+X) 



^dAu — 77~r me ja„ t-ji_/j„^/„ ; i 1/ ; v \ • (1) 

We adopt = 42 mb for the inelastic p+p cross section, while the non-elastic d+Au 
cross section, = 2.21 ± 0.09b, and the mean number of binary collisions, (A^in) = 
7.5 ±0.4, are determined by a Glauber model calculation [19]. Normalization systematic 
errors mostly cancel in the ratio. At r\ m 0, R^ Au ~ 1 , with a Cronin peak at pj > 2 GeV/c 

|0]|. In contrast, at (77) = 4.00, R$ Avi < 1, as seen in Fig. El (right). The decrease of R dAu 
with r\ is qualitatively consistent with models that suppress the nuclear gluon density 

§0,001. 

^dAu ^ S significantly smaller than a linear extrapolation of ^Au t0 ^7 = 4, 
consistent with expectations of isospin suppression of p + p — > h~ +X J4J]. 

In summary, inclusive yields of forward n° mesons from p+p collisions at \/s = 
200 GeV are consistent with NLO pQCD calculations. The cross sections show rapid 
variation with both xp and pp. The d+Au yield is consistent with a model which 
treats the Au nucleus as a CGC for forward particle production. Comparisons with 
other models will be interesting to perform. Normalizing to equal numbers of binary 
collisions, the d+Au yield is significantly smaller than the p+p yield. The r\ dependence 
of the reduction is consistent with models which suppress the gluon density in nuclei, 
in addition to exhibiting isospin effects at these kinematics. Additional measurements 
at different centralities and with other final states will help elucidate the cause of the 
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FIGURE 2. (Left) Inclusive n° production cross section for d+Au collisions, displayed as in Fig. Q 
The curves are from models which treat the Au nucleus as a CGC, normalized with a p j -independent K- 
factor. (Right) Nuclear modification factor for minimum-bias d+Au collisions versus pj. The solid circles 
are data for n° mesons. The inner error bars are statistical, the outer combine these with the point-to-point 
systematic errors, and the shaded band is the normalization error. The x's and stars are data for hr at 
smaller rj. These error bars are statistical, while the shaded boxes are the point-to-point systematic errors. 



suppression. In addition, both di-hadron correlation data and quantitative theoretical 
understanding thereof are needed to facilitate tests of a possible CGC. 
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